Feedback control of unstable cellular solidification fronts 



A. J. PonsEI and A. Karma 
Physics Department and Center for Interdisciplinary Research on Complex Systems, 
Northeastern University, Boston, Massachusetts 02115 

S. Akamatsu0 M. Newey, A. Pomerance, H. Singer|] and W. Losert 
Department of Physics, IPST and IREAP, University of Maryland College Park, USA 

(Dated: February 6, 2008) 

We present a numerical and experimental study of feedback control of unstable cellular patterns 
in directional solidification (DS). The sample, a dilute binary alloy, solidifies in a 2D geometry under 
a control scheme which applies local heating close to the cell tips which protrude ahead of the other. 
For the experiments, we use a real-time image processing algorithm to track cell tips, coupled with 
a movable laser spot array device, to heat locally. We show, numerically and experimentally, that 
spacings well below the threshold for a period-doubling instability can be stabilized. As predicted by 
the numerical calculations, cellular arrays become stable, and the spacing becomes uniform through 
feedback control which is maintained with minimal heating. 
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The control of cellular microstructures in directional 
solidification (DS) of dilute binary alloys is a subject of 
both industrial and fundamental interest Q, Q- DS is 
produced in the presence of a thermal gradient, G, which 
moves at velocity V p . Cellular microstructures arise from 
the morphological instability of a planar front when the 
velocity of the thermal gradient is above some thresh- 
old, V c , that depends on the gradient G and the alloy 
concentration. Once the planar front becomes unsta- 
ble, it restabilizes (ideally) into a periodic array of solid 
fingers or "cells" . A dynamic competition between so- 
lute diffusion, in the liquid, and capillary effects, at the 
moving solid-liquid interface, determines the typical cell 
size, but the local wavelength or cell spacing, A, admits 
a wide range of stable values. When the average cell 
spacing, Ao, is above some spacing threshold, A c , the ar- 
ray is stable and achieves a steady configuration. When 
Ao < A c , some cells, generally those of larger local wave- 
length, grow faster than their neighbors. This leads to 
the amplification of some modes and, eventually, to the 
elimination of, approximately, one cell out of two 0, Q 
(period-doubling instability). During the initial evolu- 
tion of the cell array, the elimination process is repeated, 
increasing progressively Ao, until a stable configuration 
with Ao > A c is reached. 

In a previous experimental work, Lee and Losert Q 
(see also 0) have shown, using the so called "combing 
method" , that it is possible to select a uniform cell spac- 
ing within the stable range. The combing method con- 
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sists of using strong local temperature perturbations in 
the vicinity of the front during the initial transient evo- 
lution. The thermal perturbation sets the periodicity 
to the cell array. However, a permanent control of cel- 
lular patterns outside the stability domain had not yet 
been achieved. In this letter, we propose a new scheme 
to achieve this control and demonstrate its feasibility 
in both phase- field simulations and experiments. Our 
scheme works, essentially, by slowing down the growth 
of any cell which overgrows the average position of other 
cells in the direction of growth y. For this purpose, we ap- 
ply local heating close to the protruding cell tips. A key 
feature of this scheme is that the amplitude of feedback 
perturbations (i.e. the magnitude of heating) essentially 
vanishes in the controlled state. 

We have performed numerical calculations using a 
modified version of a recently proposed quantitative 
phase-field model of binary alloy solidification 0> 0, IS 
HOj . For the experiments, we use a model transparent al- 
loy, namely, succinonitrile (SCN)-coumarin 152 (C152), 
in thin sample (see, e.g., |5j and refs. therein). 

In the numerical calculations, we use a feedback control 
scheme which is a simple step function, i.e. the amount 
of heat injected in a spot-like region is constant regard- 
less of how far the targeted cell grows beyond the average 
cell position y. We consider an array of N cells in a rigid 
box (with no-flux boundary conditions) . For a given cell 
q, local heating at the tip is applied when the distance 
from the average cell position, y q — y is larger than a 
pre-defined cut-off 8 > 0. The temperature field T(x, t), 
neglecting the production of latent heat ( "frozen temper- 
ature approximation"), can be expressed as follows: 



T{Z,t)=T + G(y-V p t)+p(Z,t), (1) 



2 



|m|coo (shift in melting temperature) 
D (diffusion coefficient) 
r (Gibbs-Thompson coefficient) 
V p (pulling speed) 
G (thermal gradient) 
do (capillary length) 
k (partition coefficient) 
£4 (0.7 % anisotropy) 



2 K 
10" 9 m 2 /s 
6.48 x fO" 8 Km 
32 iim/s 
f43.587 K/cm 
\xm 
0.3 
0.007 



TABLE I: Parameters for impure succinonitrile used in phase- 
field simulations jll|. See also |10|. 



t)=Y, ffAo H [y q (t) - y(t) - 6] exp 



-(x-x g (t)) ; 



(2) 

where x is the two-dimensional vector position, p(x, t) is 
the imposed thermal perturbation, H[...] is the Heaviside 
step function, and the sum extends to the N cells. These 
equations enter as a modification of Eqs. (132-133) in 
Ref. 10]. They approximate the fact that each heat- 
ing spot in the experiments results in the build up of a 
Gaussian bump in the thermal field of constant width, 
£, and amplitude, gA^. A precise determination of the 
phase diagram for the SCN-C152 alloy used in this study 
is presently lacking. Therefore, for sake of realism, we 
carry out the phase-field simulations for physical param- 
eters (Table |U) estimated in Ref. [ll| for a SCN-X alloy 
(where X stems for an unknown impurity). 

The feasibility of our control scheme is best illustrated 
by comparing the dynamic evolution of a strongly unsta- 
ble cellular array with and without control. The wave- 
length of the array Ao = 11.3 fxm is well below the sta- 
bility threshold A c w 50 fim for cell elimination [l^ . 
Without control, the known spatial period-doubling in- 
stability that leads to cell elimination is observed. In 
contrast, with control on, this instability is suppressed. 
This is illustrated in Figure ^ where we show the evolu- 
tion in time of the front position, y qi at specific x posi- 
tions that initially correspond to cell tips, when control 
is on. It is observed that y decreases at the same time 
that y q dispersion decreases. Later, when y q dispersion 
is reduced, the average position advances arriving to a 
steady state value. In Figure [21 we show some snapshots 
of the front evolution. Due to the large initial value of 
\Vq ~ y\ along the pattern, a relatively strong (or fre- 
quent) heating is necessary for melting protruding cells 
backwards (decreasing y). This entails a transitory disor- 
der visible in Fig. [5Jd. However, this large-perturbation 
stage ceases as soon as cell tips reach an almost equal 
undercooling (Fig. Eventually, a uniform pattern is 
stabilized (Fig. We will see below that a uniform 

A(x) distribution, obtained numerically by construction 
of the array, will also be obtained in the experiments. In- 
terestingly, the heating frequency (not shown) decreases 
during the whole process, and reaches a very small value 
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FIG. 1: Evolution in time of the tip positions, y q (t), in the 
, growth direction referred to a given isotherm with control on. 
g/G=2.5, 5/ A = 0.5% and, £/A=0.97. A=11.3 fim. 



when control has been achieved. Note also that the sys- 
tem could not be controlled by setting S — (S/A = 0.5% 
in Fig. 0). The reason for this is that it is practically im- 
possible to maintain all tips at exactly the average fixed 
position. Therefore, our feedback algorithm with 6=0 
heats continuously some of the tips and leads to a com- 
pletely melt out of cells, destroying the periodicity of the 
array. Naturally, feedback control also fails for 6 values 
larger than, say, 0.1 A. Note finally that a similar cut-off 
(about one pixel) was also introduced in our experimental 
feedback program. 

In the experiments, we solidified a binary SCN-0.1wt% 
C152 alloy in a gradient G = 10 KcmT 1 , and with 
V p = 2 — 20 /ims -1 (same device as in Ref. 5]). The 
alloy, which crystallizes in a bcc cubic crystal, is confined 
in a 100 [im thick, 2 mm wide, 150 mm long, glass-wall 
microtube. Prior to DS experiments, a single crystal with 
its [100] axis almost parallel, within less than 2°, to the 
solidification axis was selected by an empirical procedure 
and grown in such a way that it fills the container; it 
served as a permanent seed for further experiments. Ex- 
periments justify the use of a 0.7% anisotropy for the 
surface tension 13]. The local heating system is a holo- 
graphic laser tweezer system (BioRyx200 from Arryx Inc) 
which consists of a 2-W NdYAG (A=532 nm) laser that is 
focused onto a spatial light modulator. The single laser 
spot is then split into a multiple diffraction spot pattern, 
which is projected into the imaged region through a (4x) 
objective of an inverted microscope. The system software 
allows independent positioning of hundreds of spots into 
controlled arrangements. Heating in the liquid is due to 
partial absorption of light by the fluorescent dye 0] . In 
a first approximation, a laser spot (about 10 \xm in di- 
ameter) acts as a point-like source of heat which diffuses 
rapidly in a quasi-bulk medium, including the thick glass 
walls. A Gaussian bump superimposed to the linear gra- 
dient (see Eq. J2J)) is, thus, a realistic representation of 
the modified thermal field. From a rough calibration, we 
estimate that for short-exposure times (of order 1 s) and 
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FIG. 2: Snapshots of the front evolution of the system shown 
in FigureQ The near cell-tip region is displayed for simplicity. 
Times (increasing upwards): (a) 0.0 s, (b) 0.6952 s, (c) 1.3992 
s and (d) 2.1032 s. 



low (about 0.2 W) laser power, the material is heated by 
less than 0.1 Kelvin in a region which extends over less 
than 40 ^m. The numerical simulation heating power 
corresponds to <?Ao « 0.4 Kelvin and extends in £ as 10 
fim. There is no measurable overlapping between two 
neighboring laser spots, which are much smaller than Ao 
(typically 100 fxm). 

We image the solidification front with a (1024x1280 
pixel 2 ) digital camera. After some image processing, the 
front line is detected and smoothed with a moving cen- 
tral average of 5 pixels and the list of tips is detected. 
They are subsequently sent to laser controller after cal- 
culating the deviation from the average tip position. We 
use a feedback control program based on the numerical 
simulations to automatically place laser spots in the liq- 
uid region ahead of the protruding tips. By restricting 
the number of controlled spots to eight (N = 8), we are 
able to update the laser spots at approximately one Hz. 

A sequence of images during a typical feedback control 
experiment is shown in Figure [3| A simplified spatio- 
temporal diagram of the experiment is shown in Fig. 0] 
First, we perform a fast partial melting of the sample to 
obtain a planar solid- liquid interface at rest (not shown) . 
Then, we start pulling at relatively high velocity {V\ ~ 
18 /ims -1 ) which allows us to obtain a cellular array with 
small spacing by using the combing method (Fig. [3^,). 
Then, we switched off the combing laser array, turned on 
control, and decreased the velocity stepwise down to a 
final value V4 = 4.6 /ims -1 (« X^i/4) at which the pattern 
is unstable. After a relatively short transient, during 
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FIG. 3: Experimental sequence of images during feedback 
control of a small-spacing cell front pattern (time increases 
upwards). Also see Fig. 2] ( a ) Initial, high- velocity cell pat- 
tern (t = 0; V p = 18 /zms -1 ) obtained with the combing 
method (see text); (b) Transient, distorted, pattern due to 
strong perturbations just after control is turned on (t = 1100 
s; V = 6.1 /1ms" 1 ); (c) Stabilized pattern (t = 1520 s; 
V p — 4.6 /^ms _1 ); (d) Period-doubling instability after con- 
trol has been turned off (t = 1610 s); (e) Distorted pattern 
after control is turned on again (t = 1620 s); (f) Re-stabilized 
pattern (t — 2100 s); (g) Highly unstable pattern (control off) 
(t = 2300 s); (h) Large-spacing (uncontrolled) pattern. Hori- 
zontal bar: controlled region. Frames: feedback control is on. 
Horizontal dimension: 2000 /im. 



which the pattern is strongly perturbed due to frequent 
illumination (Fig. a weu controlled, small-spacing 

pattern is eventually obtained (Fig. |3J;). Note that the 
area accessible to the laser spot array is limited to about 
1 mm, thus, only one half of the cellular pattern may 
be controlled. The other half was imaged but grew freely 
without control. As it can be seen in the spatio-temporal 
diagram of Fig. many cells are eliminated outside the 
controlled window, while the number of controlled cells 
remains constant in the controlled area. 

We then switched intentionally the laser light off (af- 
ter about 10 minutes of successful control) during a few 
minutes, and observed the onset of the period-doubling 
instability (Fig. By turning feedback control on 

again, we were able to prevent cell elimination and to 
restabilize a small-spacing cell pattern similar to the ini- 
tial one (Fig. |3f). We stress the striking resemblance 
between numerical (Fig. |2J) and experimental (Figs. 01 
to [3f ) runs -including the transient stage after the laser 
is turned on (Fig. |2t> and Fig. Finally, we switched 
feedback control off again and let the instability fully de- 
velop, and, as expected, approximately one cell out of 
two is eliminated in the previously controlled area (Figs. 
and|3Ji). Let us make two important remarks. 
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FIG. 4: Lateral cell tip position i as a function of time 
(same run as Figure [JJ. Vertical bars: velocity jumps 
(Vi = 18.4 fims' 1 ; V 2 = 9.2 (Urns -1 ; U 3 = 6.1 /ims -1 ; 
Va = 4.6 /ims -1 ). Upward (downward) arrows: feedback con- 
trol on (off). Lateral drift of the pattern is due to a slight 
misalignment of the [100] axis of the single crystal with axis 

y- 

In numerical as well as in experimental runs, the spac- 
ing distribution in a well controlled pattern is remark- 
ably uniform (it is not so outside the controlled area) as 
seen, for example, in Figure [3t>- In addition, the con- 
trolled small-spacing cells sit slightly behind, thus have 
a larger undercooling than the large uncontrolled ones, 
as expected. This is not due to the added heat from the 
feedback control, but due to interactions between cells. 

The second remark concerns the heating power ap- 
plied to the cell pattern, or, equivalently, the laser spot 
exposure frequency / ahead of cells. Our main obser- 
vation is that / gradually decreases as control goes on. 
In an experiment performed at V p = 2.9 /ims -1 (not 
shown), we measured / ~ 0.5 s _1 at the beginning of 
the run, and, / « 0.05 s~ x , 10 minutes later. As ex- 
pected, the final characteristic time 1// « 20 s is smaller 
than the characteristic amplification time of the period- 
doubling instability (r w 80 s) measured in situ in an 



uncontrolled, unstable pattern. Furthermore, we always 
performed feedback control sequences much longer (more 
than 10 minutes) than r. It is interesting to note, also, 
that the final 1/f is comparable to the so-called dif- 
fusion time r d = D/V 2 = 10 - 100s (D falls in the 
10~ 10 — 10 -9 m 2 s~ 1 range), which signals a quasi steady 
regime. Therefore, control can be maintained as long as 
wanted. 

In conclusion, we have shown experimentally and nu- 
merically that highly unstable cellular solidification ar- 
rays can be stabilized using feedback control. Experi- 
ment and simulation showed similar dynamical evolution 
as the array was controlled. Unlike feedback stabiliza- 
tion schemes for planar fronts slightly above V c proposed 
by Savina et al, |15| . our approach stabilizes efficiently 
highly unstable states and can be implemented exper- 
imentally with only a discrete number of controllable 
heating points. An alternate control scheme which in- 
cludes a tunable heating pow er for each spot has been 
implemented numerically jlfij : an experimental realiza- 
tion is currently in progress. Although a quantitative 
comparison between numerical and experimental results 
is not yet possible because of the lack of detailed knowl- 
edge of the alloy phase diagram, the observed qualitative 
behaviors are identical. 

The use of other localized physical perturbations gen- 
erated e.g. by x-rays or ultrasound could potentially 
make it possible to extend this control scheme to metal- 
lic alloys where the growth of much finer array struc- 
tures with superior mechanical properties is of consider- 
able practical interest. In addition, the ability to access 
experimentally unstable steady-state patterns provides 
an important new tool to enrich our fundamental under- 
standing of pattern formation in directional solidification. 
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